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Abstract. We prove a global foliation result, using areal time, for T 2 symmetric spacetimes with 
a positive cosmological constant. We then find a class of solutions that exhibit AVTD behavior 
near the singularity. 
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. 1. Introduction 

<D ■ 

fX^ ■ Strong Cosmic Censorship (SCC), one of the outstanding conjectures of mathematical relativity, 

claims that a generic maximal globally hyperbolic solution of the Einstein equations (vacuum, or 
for a specified matter model) cannot be extended across a Cauchy horizon. Model studies of SCC 
in certain families of spacetimes with specified isometries — e.g., the spatially homogeneous solutions 
[TO] , and the T 3 Gowdy spacetimes [IS] — show that, as a step toward proving SCC in such families, 
it is very useful to investigate whether the solutions exhibit either local oscillatory ( "mixmaster" ) 
behavior near the singularity, or the more special asymptotically velocity-term dominated ( "AVTD" ) 
behavior in that region. These studies also show that the existence of a geometrically motivated co- 
ordinate system compatible with the singularity, such as the "areal coordinates" of the T 2 -symmetric 
spacetimes [3], is very useful. In this paper, we show that T 3 Gowdy and T 2 -symmetric solutions 
with a positive cosmological constant admit areal coordinates covering their globally hyperbolic re- 
gions. We also use Fuchsian methods to show that there are sets of spacetimes of this sort with 
AVTD behavior. 

Since the point of this work is to show that a number of the results which have been proven for 
T 3 Gowdy and T 2 -symmetric vacuum solutions also hold when a positive cosmological constant is 
added, we first recall the results which hold for these two classes of vacuum spacetimes. For the 
T 3 Gowdy solutions, which are defined by the presence of a U(l) x U(l) (or T 2 ) isometry group 
acting on the T 3 spacelike foliations, and by the vanishing of the twist constants, the following has 
been shown: (i) Global areal coordinates (for which the constant time hypersurfaces are defined to 
be those spacelike hypersurfaces which contain T 2 orbits of equal area) exist for all solutions [13] . 
with the time coordinate covering the range (0, oo) except in certain special (Kasner) cases, (ii) 
Generically, the solutions exhibit AVTD behavior [17J . (iii) Generic solutions have their curvature 
unbounded as one approaches the singularity, and cannot be extended outside the maximal globally 
hyperbolic spacetime region [17]. Hence SCC has been proven for the T 3 Gowdy solutions. 

For the T 2 -symmetric vacuum solutions with non vanishing twist constants, less is known. How- 
ever, it has been shown that all such (maximal, globally hyperbolic) solutions admit global areal 
coordinates [3], again with the time covering the interval (0, oo) except for certain special (Kasner) 
cases [TT]. It has also been shown [3] that, if one restricts one's attention to polarized T 2 -symmetric 
vacuum solutions, in which one of the gravitational degrees has been turned off, there is a large col- 
lection of such solutions which show AVTD behavior. It has not yet been shown that these polarized 
solutions generically exhibit AVTD behavior; however based on numerical and heuristic studies, it 
is expected that this is the case. Similar such studies [4J suggest that in the non polarized case, the 
solutions generically show mixmaster behavior. 

What happens to these results if one considers spacetimes which satisfy the Einstein equations 
with source fields? For the special case in which the source fields are collisionless (Vlasov) matter 
fields, much is known: It has been shown that both the Gowdy [1] and the T 2 symmetric [2] solutions 
with Vlasov matter admit areal coordinates, and strong cosmic censorship has been proven [7] in 
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both cases. These results depend on the special features of Vlasov matter, and the techniques are 
quite different from those used in the vacuum case. 

Motivated by recent observations suggesting that a large fraction of our universe consists of a 
component with negative pressure [19] , and noting that a positive cosmological constant is one pos- 
sible explanation of this |15j , we begin in this paper to study what happens if a positive cosmological 
constant is present (with no Vlasov fields). We obtain results similar to those for the T 2 symmetric 
vacuum solutions. In particular, we show that T 2 -symmetric spacetimes with a positive cosmologi- 
cal constant admit an areal foliation, with the time coordinate taking values in (c, oo), where c is a 
nonnegative number. In addition, we find a class of T 2 -symmetric solutions with positive cosmolog- 
ical constant that exhibit AVTD behavior near the singularity. We note that the spacetimes with 
AVTD behavior that we discuss here are not fully polarized, but also have less than the maximum 
number of free functions. Analogous to the [/(l)-symmetric solutions with AVTD behavior studied 
in [10] , they might be described as "half polarized" . The extension of AVTD results from polarized 
T 2 -symmetric solutions to half-polarized T 2 -symmetric solutions with cosmological constant which 
we show here does not in fact rely on the cosmological constant being non zero; so our results here 
strengthen the vacuum results of [S]. 

The proof of the general, global existence of areal coordinates was first presented in [5] for the 
vacuum T 2 -symmetric case. In our discussion here, we follow the slightly different approach of pQ, in 
which the global areal foliation is proved for Gowdy spacetimes with Vlasov matter. It is reasonable 
that the Vlasov case would serve as the most convenient model, since adding a cosmological constant 
is equivalent to a particular choice for the stress energy tensor. We do not here describe the complete 
proof, but rather we focus on the parts of the proof which are specifically affected by the presence 
of the cosmological constant. 

To show AVTD behavior we follow Rendall and Kichenassamy in utilizing the Fuchsian algo- 
rithm. An introduction to this method can be found in [12j and [9] . In Section 4 we briefly review 
the algorithm and identify a set of analytic initial data whose development has AVTD behavior near 
the singularity. We note that it follows from the work of Dafermos and Rendall [5] that T 2 sym- 
metric spacetimes-vacuum, with Vlasov matter, or with positive cosmological constant-are always 
inextcndiblc in the expanding direction. Hence to prove strong cosmic censorship for T 2 symmetric 
spacetimes, it is sufficient to understand the behavior of the solutions near the singularity. 

2. Field Equations 

The T 2 -symmetric spacetimes, like the Gowdy spacetimes, admit a T 2 isometry with spacelike 
orbits; but unlike the Gowdy spacetimes, they have nonvanishing twists. Letting X and Y label the 
Killing vector fields that generate the isometry, we find that the twists are defined as 

K x = e a pisX a Y f) V y X s and K Y = e a p lS X a Y p V<Y s . (2.1) 

As in [3] we rotate X and Y so that Kx is zero and Ky = K is not. This can be done without loss 
of generality. The Cauchy surfaces of these spacetimes (with non vanishing twist) necessarily have 
T 3 topologjQ, with the coordinates (x, y, 9) all having period one. The isometry is in the x and y 
directions where the above rotation has already been done, leaving 9 as the only space direction in 
which the fields can vary. Thus, the isometry reduces the problem to a 1 + 1 PDE system, greatly 
simplifying the analysis. 

It is convenient in the global areal foliations proof to use two different coordinate systems. 
Conformal coordinates are used to the past of the initial surface (towards the singularity) and areal 
coordinates are used to the future of the initial surface (in the direction of expansion). In conformal 
coordinates the area of the group orbits, R, is a function of 9 £ S 1 and (eK. In areal coordinates 
we set the area function R — t; hence the name "areal" coordinates. The field equations we consider 
are the Einstein equations with a positive cosmological constant A, 

G[a> + A 5m „ = 0. (2.2) 



lr The Gowdy spacetimes, with T 2 isometry and vanishing twist, can have Cauchy surfaces with T 3 , S 3 , or S 2 X S 1 
topology 8 . Those Gowdy spacetimes discussed here, as well as in most discussions of strong cosmic censorship in 
Gowdy solutions, have the T 3 spatial topology. 
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In areal coordinates the metric takes the form 

g = e 2 ^- u) (-adt 2 + d9 2 ) + e 2U [dx + Ady + (G + AH)d0} 2 

+ e- 2U t 2 [dy + Hdd} 2 . (2.3) 

Here U, A, a, i], G, and H are functions of 9 £ S 1 and t £ M. + . Note that the functions labeled Mi 
and M2 in [3] have been set to zero by an appropriate coordinate change. The Einstein equations 
then are as follows 

Constraint equations: 

7 } = U? + all 2 + ^(A 2 + aA 2 ) + ae 2 ^A + ^aK 2 , (2.4) 
m 

*l = WtU , + - A ,A,--±, (2 .5) 

a t = -^aV^-^A - ~pO?K 2 . (2.6) 
Evolution equations: 



a g r)e a t r) t a 2 g a ee „ 2 . Tj2 
+ e —(A 2 -aA 2 )+aAe 2 ^-^K 2 



(2.7) 



Uu - a U 0e = -^ + ^ + ^+ e w {A 2 - aA 2 ) + oAeXHO, (2.8) 

A H - aA ge = + + ^ - iAtUt + AaAeUg. (2.9) 

t 2 2a 

Auxiliary equations: 

= G t +AH t , (2.10) 
e 2 '' 

= H t - Va~^-K. (2.11) 

In conformal coordinates , where R is a function of 6* and t, we can set a(9,t) = 1. The 
metric then takes the form 

.g = e 2 ^^ (-dt 2 + d# 2 ) + e 2U [dx + Ady + (G + AH)d9] 2 



4U 



e 



- 2U R 2 [dy + Hd9} 2 . (2.12) 



The Einstein equations take the form 
Constraint equations: 

= U 2 + U 2 + ^(A 2 + A 2 ) + - ^ - ^ + Ae 2 ^) + |^ 2 , (2.13) 

= 2 W + ^A t A e + - ^ - 2^. (2.14) 
Evolution equations: 

Cfe - Cto = ^ - ^ + fJ(A 2 - A3) + Ae 2 (^>, (2.15) 
A tt -A ee = ^-^ + 4(A e C/ e _A t [/ t ), (2.16) 

P 2»J 

i? tt - = 2i?Ae 2 ("- l/ ) + ^ if 2 , (2.17) 



4U 



3e 2 " 



- m = L/ 2 - £/* + ^(A 2 - A 2 ) + A*«»-W - ^ K 2 . (2.18) 
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Auxiliary equations: 



Q = Gt + AH t , 

e 2 " 



(2.19) 




(2.20) 



We recover the equations for the vacuum Gowdy spacetimes by setting the cosmological constant 
and the remaining twist constant K to zero, and letting both a{6, t) — 1 and R(6,t) = t. In this 
situation the constraints decouple from the evolution equations and can be effectively ignored in the 
course of the analysis until the end. Adding any one of the following — nonzero twist, cosmological 
constant, or Vlasov matter — does not allow us to have both a = 1 and R = t. Hence, the constraint 
and evolution equations are necessarily coupled in this setting. 

Local existence for each of these coordinate systems is proven in [5] . 



Let 7 be a smooth Riemannian metric on T 3 and let tt be a smooth symmetric tensor on T 3 which 
together satisfy the Einstein constraint equations with positive cosmological constant A and which 
are invariant under an effective T 2 group action. We then call (T 3 , A; 7, 7r)-which we shall sometimes 
abbreviate as (7, ir)-& set of T 2 -symmetric Einstein-A initial data . The following theorem, in the 
spirit of [3] and p] , is our main result on the existence of global areal coordinates for the spacetime 
development of such data: 

Theorem 1. Let (T 3 , A; 7,71") be a set of T 2 - symmetric Einstein-A initial data on T 3 . For some 
nonnegative constant c, there exists a globally hyperbolic spacetime (M 4 ,<?) such that 



(2) g satisfies the Einstein equations with positive cosmological constant A. 

(3) M 4 is covered by areal coordinates (x,y,8,t) with t € (c, 00), so the metric globally takes the 
form ([23]) . 

(4) (M A ,g) is isometrically diffeomorphic to the maximal globally hyperbolic development of the 
initial data (T 2 , A; 7, 7r). 

Outline of the proof. The basic steps of the proof follow the pattern established in {3] and continued 
in [TJ. In particular, the first part of the analysis involves the study of the fields to the past of the 
initial surface, in the contracting direction. Working in terms of conformal coordinates, we show 
in Section 3.1 that as long as the area function R(8,t) stays bounded away from zero, the past 
maximal globally hyperbolic development in terms of conformal coordinates for this initial data, 
denoted £'^ )n f(7, tt), has t — > —00. 

The next part of the proof involves the establishment of a number of geometric results related to 
the behavior of the function R in the spacetime development of the initial data. One shows that R is 
positive everywhere in the maximal development D{pj, n) of the initial data, and that R approaches 
a nonnegative value Rq along any past-directed inextendiblc timclikc path contained in D~ J^, tt). 
Further, one finds that this limit Rq is the same limit along all such paths. Thus one can identify Rq 
with the number c which appears in the statement of Theorem 1. One now proceeds to show that 
^conf(7' 7r ) admits a foliation by constant R Cauchy surfaces, and in addition one shows that if the 
conformal time t — > —00, then D~ J^, ir) is identical to the past maximal development D~(~/, n) of 
the initial data. All of these geometric results are proven in Section 5.2 of [3] for the vacuum case. 
The proofs are essentially unchanged by the addition of the cosmological constant, so we leave them 
out of this paper, and refer to [3|. 

The proof in the expanding direction is more direct, as we work in areal coordinates. We choose 
one of the R = constant Cauchy surfaces from the past of the initial surface. The spacetime induces 
data (71, 7Ti) on this surface. The future maximal globally hyperbolic development in terms of areal 
coordinates for this data, denoted-D^. eal (7i,7Ti), has t — > +00. This is shown here in Section 3.2. To 
complete the proof, we cite another result from [3] which shows that -D^. ea i(7i, tti) ps D + (71,711). 



3. Global Areal Foliation 



(1) 



M 4 = T 3 



x (c, 00) 
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3.1. Contracting Direction. In this section we sketch the steps of the first part of the proof, 
which establishes global existence in terms of conformal coordinates in the contracting direction, 
so long as the metric function R stays bounded away from zero. We follow the overall outline of 
the proof carried out in pQ, detailing those steps in which the presence of a non zero cosmological 
constant A plays a significant role. 

We recall that the local existence of a solution in conformal coordinates has been demonstrated 
by Chrusciel in [5]. To show that the past development of the initial data extends to t ~ * — oo so 
long as R stays bounded away from zero, it is sufficient (as a consequence of Theorems 2.1 and 2.2 in 
[13j ) that we obtain uniform C 2 bounds on the metric functions. To verify these bounds, we proceed 
as follows: 



Step 1. (Monotonicity of R and C 1 bounds on R.) 

We first note that we can use the arguments of Theorem 4.1 of [5] to show that Vi? is timelike, 
i.e., that g(\7R, \7R) < 0. Those arguments, which in [5] are carried out for T 2 -symmetric solutions 
with A = and non zero twist, apply here with A > since they rely essentially on the positivity 
of the twist term in the A = version of the constraint (2.13), and not on the explicit form of this 
term; hence we find that the A term in (2.13) has the same effect. An immediate consequence of 
this is that Rt must be nonzero everywhere. We choose Rt > so that the past does in fact coincide 
with the contracting direction. 

Next we define the null vectors 

d x = -j=(d t -d 8 ); dt = -j=(dt + de), (3.1) 



and from the evolution equation for R we calculate 

n ± 

4i? 3 

or, equivalently 



R H - RAe 2 ^ + ^K 2 , (3.2) 



P 2v 

Rtx = RAe 2 ^'^ + -^K 2 . (3.3) 

Using A > 0, we see that the right hand side of equation (|3.2|) is positive. We can then integrate 
back along the null directions, A and £, starting at (0q, to) on the initial surface, and verify that 

Rt(0 o + s,to-s) <Rs(0 ,t ), (3.4) 
Rx(0o-s,t o -s)<R x (9 o ,t o ). (3.5) 

Thus, for any t < to and 8 £ S 1 we have 

R ( (t,0) < sup R s (t Qr ), (3.6) 

6 

Rx(t, 6) < sup Rx(t ,-). (3.7) 

6 

Since R t = 1/V^(R^ + R\), it follows that 

Rt(P,t)<^j=8uv(Rt + Rx)(to,-)- (3-8) 
V2 e 

Hence we have Rt bounded to the past. Since Vi? is timelike, \Rt\ > \Rg\, so we have have established 
C 1 bounds on R. Noting that in deriving these results we only use the the sign of R^x and Rx^, it 
follows that the addition of a positive cosmological constant does not change the argument used in 
the vacuum T 2 -symmetric case. 



Step 2. (Bounds on U, A, and their first derivatives.) 

We define here the same energy-like quadratic functions G and H which are used in the vacuum 
T 2 -symmetric case [3] and in the Gowdy with Vlasov matter case [T] . (There is no relation between 
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these forms and the metric functions of the same names.) 



' u 2 ) ' 



AU 



G - -R{Uf + Ui) + —(Af + Aj), (3.9) 



AR 

Taking derivatives along the null directions we find 



H = RUtUe + ^—AtAe. (3.10) 



l p 4U 

d x (G + H) = RUsAe^-u) _ _ ^2 _ ^ + —(-A* + Aj)], (3.11) 

1 p 4U 

c%(G -H) = RU x Ae 2 ^ - -R X [U? - U 2 + ^(-A* + Aj)]. (3.12) 

If we now integrate these equations along null paths, starting at (t, 9) with t < and ending on 
the io-surface, we obtain 

[G + H](t ,e- (t -t)) - [G + H](t,6) = £ ° K x + U s Td\, (3.13) 

[G -H}(t ,e+ {t - t)) - [G - H]{t, 9) = j ° K 2 + U x Td£, (3.14) 

where 

1 e 4U 

Ki = - U 2 e + ^ (~ A t + A e)l (3-15) 

1 p 4E/ 

K 2 = -^Rx[U 2 - U 2 + ^(-^ + A 2 e % (3.16) 

T = RAe 2 ^- u \ (3.17) 

To handle the T terms in (|3.13[) and (|3.14[) . we proceed as in the Vlasov case [T]: Recall from 
the evolution equation for R that we have 

R H = RAe 2 ^ + f^K 2 , (3.18) 
and a similar equation for R^ x . Integrating these along null paths we find that 

/to „2r; 

[RAe^v-u) + —K 2 ](s, 9+(s- t)) ds, (3.19) 

/to „2r; 

lRAe 2 ^-u) + -^K 2 ](s, 9-(s- t)) ds. (3.20) 

Since we have C 1 bounds on R (from Step 1), the left hand side of each of these equations is bounded. 
Hence the right hand side is also bounded. This implies further that the term 

rto rto 



J ° [T] (s) ds = [RAe 2{, '- u) }(s) ds (3.21) 



is also bounded on the interval (t,to] x S 1 . 

To handle the integrals of the K\ and K 2 terms in (|3.13[) and (|3. 14[) . we proceed as in the 
vacuum case [3] by bounding them in terms of the quadratic form G in the following manner. From 
the definitions of K\, K 2 , and G as well as from the C 1 bounds on R, we have 

■ , _ . CG , § § CG 
\Ki\<— and \K 2 \<—. 

Then using the elementary inequalities ±2ab < a 2 + b 2 , we also have 

|C/d<(^) V2 and \U X \<(^)V\ 

We now add equations (|3.13[) and (|3.14[) and take the suprema over 9. Using various standard 
inequalities along with the crude estimate \J~G < 1 + G, we obtain an inequality which permits the 
application of Gronwall's Lemma. This gives us a uniform bound on sup e G(t±, •), provided R stays 
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bounded away from zero. This bound on G(t, 9) gives us bounds of \Ut \ and \U$\, which we can then 
integrate to get bounds on U itself. We then have bounds on \A t \ and \Ag\, which can themselves be 
integrated to obtain a uniform bound on A. Hence we have C 1 bounds on U and A. We note that 
the terms introduced by the presence of the (positive) cosmological constant are easily controlled by 
our C 1 bounds on R. 

Step 3. (C 1 bounds on rj.) 

As in [3] we define the quantity 

p:=r)+^ln(R). (3.22) 

Then using the evolution equations for r) and R we derive 

/3 A? = Z + 2Ae 2{ri ~ u \ (3.23) 

where Z contains only quantities which we have previously shown to be bounded. Integrating this 
along null paths, along with the corresponding equation for (5^\, we find that 

Px(t o ,O + (t o -t))-p x (t,0) = f [Z + 2Ae 2 ^-^}( Si e+( S ~t))d S , (3.24) 



together with a similar equation for (3^(t,9). We now recall from an argument presented in Step 2 
that, given C 1 bounds on R, we determine that the integral 



J RAe 2 ^-^ dl; 



is bounded. It follows that (3\ is bounded on (t, to] x S 1 provided R > e > 0. The same argument 
applies to Since (3 t — 1/v / 2(/?a + /%) and (3g — 1/v2(/3a — both (3 t and (3g are bounded on 
(t, to] x S 1 . Given the definition of f3 and the C 1 bounds on R, we see that rjt and rje are bounded 
on (t,to] x S 1 . Integration gives a bound on rj itself. 

We now have C 1 bounds on the principle metric functions. These bounds on U and r], together 
with the condition that R is bounded away from zero, allow us to control all terms which contain a 
cosmological constant, since they always appear in the form 

Ae 2 ^) or RAe 2 ^- u \ (3.25) 

Thus the remainder of the proof in the contracting direction is unaffected by the presence of the 
cosmological constant and follows as in [3]. 

3.2. Expanding Direction. The work we have just described shows that, for the initial surface 
and T 2 -symmetric data we have chosen on it, the past (contracting direction) of this initial surface 
is foliated by R — constant slices. Choosing one of these slices as the new initial hypersurface (and 
labeling its R = constant time coordinate as R — t^ we now have induced initial data appropriate for 
areal coordinates, and we seek to show global existence of the evolution in terms of these coordinates 
in the expanding direction, to the future of ti- As in the argument for the contracting direction, 
the key to showing global existence is to establish C 2 bounds for solutions of the areal form of the 
equations (|2.4|) - (|2.9p . We discuss how to do this in the following steps, focussing on those steps in 
which the presence of a positive cosmological constant plays a role. 

Stepl. (Bound on \U\.) 

We first establish monotonicity for the energy function E(t), defined by 



E(t) 

where 



si 



At 4 



d9, (3.26) 



IT 2 p 4U A 2 

h= t + ^ u 2 + { t +V ^ A 2y (3 27) 

Using the evolution equations for U and A and integrating by parts we determine that 



s 
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We note that if we leave the K 2 term out of the definition of E(t), then this new function is also 
positive and monotone decreasing. However, if we omit the A term, then we no longer have a 
monotone function, as a consequence of the presence of A in the evolution equation for U. The 
situation is analogous to the Gowdy-Vlasov case, in which the mass energy density quantity p 
appears in the definition of the energy function E. 

We now use this behavior of E(t) to establish bounds on \rj\ and \U\. Defining the new function 
/3by 

0:=r) + ±1n(a), (3.29) 
we use the constraint equations for r\ and a to show that 

fit = t[U 2 + all 2 + y(A* + aA 2 )} - taAe 2 ^ - e —^-. (3.30) 
We then immediately see that 



p t <t[U 2 + aU 2 + e w {A 2 +aA 2 )], (3.31) 



4(7 



e 2 ^aK 2 




t > -taAe 2 ^ - ^-pt-. (3.32) 

Clearly the presence of the positive cosmological constant does not affect the sign of the second 
inequality, (j3.32|) . The monotonicity of E(t), together with equation (|3 . 3 1 1) . is used to get an upper 
bound on (3. We refer the reader to [3] for the details, as the argument is no different from the 
vacuum case. (Note, however, that v — —0 is used in [3] so the argument for an upper bound of 
here follows that for a lower bound of v in [3].) Unlike the vacuum case, however, to verify a lower 
bound on we must first bound \U\, as can been seen from the A term in (|3. 321) . We accomplish 
this as follows: (See also [1]). From Holder's inequality we have 

\U(t,9 2 )-U(t,6 1 )\ = | / ~U e {t,9)d0\ 

\ 1/2 /A _ \ 1/2 

1/2 dd ) [J y/aUeffl) ■ ( 3 - 33 ) 

The second integral in (|3.33p is clearly bounded by (Efo)) 1 / 2 - To bound the first integral we use 
the evolution equation for a to find 

dtia- 1 ' 2 ) = 2t^e 2 ^-u) A + e2 "^ g2 , (3.34) 

Integrating this over time and space we have 

f\- 1 ' 2 {t,6)de= f[2s f 2 yfce^-^k + e2 "^ 2 ] d6ds + /' a^ih, 6) dQ 

J6i Jti JOi ^ s JSi 

< [ 2sE(s) ds + 2-kC 

< E(t 2 ){t 2 -t\) + 2irC. 
Thus, for all t £ ^2^3) an d 9\, 62 € S 1 , we have 

\U(t,e 2 )-U(t,9 1 )\ <C(t), (3.35) 

for some bounded function C(t) of t. This controls the oscillation of U(t,8). To control \U(t,9)\, we 
proceed as follows. Note that in the following, C(t) changes from one line to the next but always 
refers to a bounded function. 
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We now consider 

U(t,9)d6\ = \ f f U t (s,8)d8ds + f U(t 2 ,6)d6\ 
s 1 Jt 2 Js 1 Js 1 



< f f \U t (s,8)\d8ds + \C\ 
Jto Js 1 



l[i x>><J9) / ( I ^jLde) ds+\C\. (3.36) 



'a 



Again, the second integral in (|3.36[) is bounded by (Efa)) 1 / 2 ■ Also, as a consequence of the constraint 
equation for at, a is bounded and hence y/a is bounded. As before, the monotonicity of a depends 
on our assumption of positive A. Thus we have the inequality 

U(t,8)d9\ <C{t), (3.37) 

for some bounded function C{t). We now use (|3.35[) and (|3.37[) to get a bound on \U\. If we define 
U + := maxfl U{t, •), we find 

U+(t)= [ U{t,8)d6+ [ (U + (t)-U(t,6))d6. (3.38) 
Js 1 Js 1 

Using (|3.35p . where 8 2 is the value for which U achieves its maximum and 9± = 6, we find 

\U + (t)-U(t,8)\<C(t) 
for all 8 £ S 1 . Taking the absolute value of (|3.38[) we find 

\U+(t)\<\[ U(t,8)d8\+ [ C(t)d8. (3.39) 
Js 1 Js 1 

Using (|3.37p . we have bounded |C/+(£)| on [£2, £3)- Repeating this procedure for U-(t) :— mine U(t, •) 
we arrive at a bound on |{7_(t)|. Thus, \U\ is uniformly bounded on S 1 x [£2,^3)- 

Step 2. (Bounds on /3, a, a t , and \rj\.) 

Given the upper bound on (3 and the bound on \U\, we see from equation (|3 . 32[) that 

- C(t) < (i t (3.40) 

for some bounded function C{t). Integration gives 

t 

-C(s)ds + f3(t 2 ,8) < (3(t,8), (3.41) 

and f3(t,8) is bounded from below. Thus \f3\ is bounded on S 1 x [t2:t^). 

Bounds on a, at, and r\ follow almost immediately. Recall the constraint equation for at, written 

as 

dt i n (a) = 2* = -4tAe 2 ^ - ^£ (3.42) 

The bounds on |/3| and |J7|, after an integration, give us a bound on | ln(a)|. Hence, a is bounded 
away from zero. The constraint equation for a t immediately gives us a bound on a t . Finally, from 
the definition of j3, 

/3 = 7?+i In (a), (3.43) 

is is clear that |??| is also bounded on S 1 x [£2, £3)- 

The metric functions \U\, and |?7| are bounded on S 1 x [i 2 , £3) and a is bounded away from zero 
there. In areal coordinates, all terms induced by the cosmological constant have the form 

ae 2(v-u) A Qr 4ta 2 e 2( "^A. (3.44) 

Whenever such terms arise in the remainder of the proof they are bounded and the argument 
proceeds as in [3]. 

This completes our proof of Theorem 1. 
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4. Analysis in the Neighborhood of the Initial Singularity. 

One of the prime motivations for showing that a spacetime admits areal coordinates globally is 
so that one can use the coordinates to locate where singular behavior occurs. In particular, if one 
can show that a spacetime admts areal coordinates with t running from zero to infinity, then one 
knows that there is a big bang type singularity (whatever its nature) at t — > 0. One can then study 
its properties by studying the metric, in the form (|2.3[) . in the neighborhood of t — > Q. 

We have just proven that every T 2 -symmetric solution with positive cosmological constant admits 
global areal coordinates, but have not shown that the t coordinate extends to zero. Now in the 
vacuum T 2 -symmetric case (with A = 0), it has been proven [llj that t covers (0,oo) in almost all 
cases. The one exception is when the initial data is flat Kasner; then one may have t > c > 0. In the 
T 2 -symmetric case with Vlasov matter [5D] , and in the vacuum Gowdy case , we also know that 
the areal coordinate t covers all of (0,oo). While we have not yet determined whether the addition 
of the cosmological constant changes this, all evidence suggests that it does not, and so we shall 
presume in what follows that we are working with spacetimes in which t — > 0. 

We wish to show here that there are T 2 -symmetric solutions with positive cosmological constant 
which exhibit AVTD behavior near the singular region at t = 0. In the course of doing this, we shall 
in fact find that there is a much wider class of such spacetimes than expected; further this widening 
of the class of AVTD solutions holds for T 2 -symmetric solutions whether or not the cosmological 
constant is nonzero. Thus we extend the earlier results obtained in [9]. 

To do this, we use the Fuchsian algorithm, which has in previous work been successfully employed 
to show that there are AVTD solutions in the vacuum T 3 Gowdy family and in the polarized T 2 - 
symmetric family of solutions. To understand why the Fuchsian method is effective, we recall what 
AVTD behavior is: Roughly speaking, this condition holds if, near the singularity, terms involving 
spatial derivatives of the fields become less important in the evolution equations than those involving 
time derivatives. Thus, each spacetime observer approaching the singularity effectively "decouples" 
from other observers and evolves in their own "private universe" near the singularity. In many cases 
these private universes are locally Kasner spacetimes with each observer having their own Kasner 
exponents. 

The Fuchsian algorithm is well adapted to verifying AVTD behavior in cosmological spacetimes, 
because it is designed for analyzing solutions of certain types of differential equations in the neigh- 
borhood of singularities, and more specifically for showing that these solutions are bounded and 
well-behaved as one approaches the singularity. Hence one may proceed as follows: One first forms 
a new system of PDEs (labelled here the "truncated system" ) , which is obtained by eliminating 
the spatial derivative terms from the areal form of all of the Einstein equations except for those 
which contain nothing but spatial derivative terms. Next, one finds functions (a parametrized set 
of them) which are asymptotic solutions to the truncated system, in the sense that as t — > these 
functions approach solutions to the truncated system (these functions are called the "asymptotic 
data" ) . One then writes the metric fields as sums consisting of the asymptotic data plus remainder 
terms, and one substitutes these sums into the complete Einstein system. The resulting equations 
act as a parametrized system of PDEs for the remainder functions, with coefficients depending on 
the asymptotic data. If these equations can be written in the "Fuchsian form" 

(td t + M(9))v = t e f(t, 9, v, d e v), (4.1) 

where the vector v represents the remainder terms, where M{9) is a matrix which depends on the 
spatial coordinate 9 and satisfies certain positivity conditions (noted below in Step 3), where e > 0, 
and where / satisfies certain analyticity conditions (also noted in Step 3 below), then a solution v 
exists near t —> zero, with v vanishing as t approaches 0. It follows (if indeed we do achieve Fuchsian 
form) that there are solutions of the full system of equations which approach solutions of equations 
in which spatial derivatives have been neglected. We thus have solutions of the Einstein equations 
with AVTD behavior, parametrized by the choice of the asymptotic data. 

Our experience with vacuum T 2 -symmetric spacetimes strongly suggests that they do not gen- 
erally exhibit AVTD behavior; but rather show the more complicated Mixmaster behavior (see [4]). 
The presence of a positive cosmological constant is not expected to change that. Our previous work 
with the vacuum case [9] suggests that AVTD behavior is to be found in polarized T 2 -symmetric 
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solutions, which are characterized by setting two of the initial data functions to constants. In this 
work, besides extending such results to the A > case, we also show that indeed the class of solu- 
tions with AVTD behavior is more general than the polarized solutions: Both with and without the 
cosmological constant, we obtain (using Fuchsian techniques) solutions with AVTD behavior which 
are "half-polarized" in the sense that only one of the asymptotic initial data functions is set to a 
constant. This condition is stated more specifically below. 

To set up our results, looking toward obtaining equations in Fuchsian form, it is useful to first 
rewrite the areal coordinate version of the Einstein equations (2.4)-(2.9) using the the operator 
D := td t - (We do not include equations (2.10)-(2.11), since these are essentially decoupled from 
(2.4)-(2.9), and thus can be handled at the end.) We obtain 



1 t 2 p au 

D 2 U - t 2 aU g = — (Da)(DU) + -a g U e + {(DA) 2 - t 2 aA 2 ) + t 2 ake 2 ^- u \ (4.2) 

1 t 2 

D 2 A - t 2 aA ee = 2{DA) + — (Da)(DA) + -a e A B - 4(DA)(DU) + 4at 2 A e U e , (4.3) 

Da = -At 2 a 2 e 2 ^-U)k - ^ 2 K 2 , (4.4) 

D V = (DU) 2 + t 2 aU 2 + ^ ((DA) 2 + t 2 aA 2 ) + t 2 ae 2 ^ A + ^K 2 , (4.5) 

m = 2(DU)Ue + e ^(DA)A e (4.6) 

Since equations (4.4) and (4.5) contain time derivatives (even though they are constraint equa- 
tions), we shall treat them as evolution equations in our analysis here. We note that the vanishing 
of the constraint function 

N = Ve -2(DU)Ue-^(DA)A e + ^. (4.7) 

is preserved in time, since it follows from (4.2)-(4.5) together with the identity d t rjg — dgrj t = that 
N satisfies the linear, first-order, ordinary differential equation 

d t N-^N = 0; (4.8) 

if N(to) = for any time to, then N(t) = for all time. 

We now carry out the AVTD argument in a series of steps: 



Step 1. (Asymptotic solutions to the truncated system.) 

Eliminating all terms with spatial derivatives in the evolution equations (4.2)-(4.5), we obtain 
the truncated system 

D 2 U = ^-(Da)(DU) + ^(DA) 2 + t 2 ake 2 ^- u \ (4.9) 

D 2 A = 2(DA) + —(Da)(DA) - A(DA)(DU), (4.10) 
2a 

Da = -At 2 a 2 e 2 ^A - ^a 2 K 2 , (4.11) 

D V = (DU) 2 + ^2 (DA) 2 + t 2 ae 2 ^- u ^A + ^K 2 . (4.12) 
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We wish to show that for a class of restricted choices of the spatial functions k(8), <f>(6), ip{9), B(9), 
X(8) and 7(0) (these are effectively functions of integration), the following 

U(0,t) = ±{l-k(9))\nt + <t>{6), (4.13) 
A(6, t) = V(0) + t 2k< - e) B(9), (4.14) 

V(8,t) = \(l-k(9)) 2 \nt + X(9), (4.15) 

a(9,t)= 1 {9), (4.16) 

are asymptotic solutions of the truncated system (|4.9p ~ (14.12p . First, substituting these choices into 
(|4.9p . we have 

= 2e^k 2 B 2 t 2k + 7 e 2 ( A -«At^ fe2 + 3 ), (4.17) 

which is satisfied in the limit of t — > 0, provided k(9) > 0. We note that it is not the A term which 
places this restriction on k(9). From equation (|4.10|) . we obtain 

Ak 2 t 2k B = Ak 2 t 2k B, (4.18) 

which holds as t — > for all values of k. Equation (|4.1ip with the functional forms (|4.13p - (|4.1Gp 
substituted in produces 

= -4Aa 2 t 1 / 2 ( fc2+3 ) e 2 ( A -^ - a 2^2 e 2A t i/2( fe -3)(fc+i) i (4 1Q) 

which is satisfied as t — > 0, provided k(9) > 3. Again, it is not the A term which places the restriction 
on k, but rather the twist term K 2 . Finally from (|4.12p we derive 

1 (1 _ k f = 1 (1 _ k) 2 + f,2 B 2 e ^ t 2k + Aat l/2(fe 2 +3) + ^! e 2A t l/2( fc -3)(fe+l) 5 (4 2Q) 

which is also satisfied for k{9) > 3. 

We conclude that the parametrized functions (|4. 13|) - (|4. 16[) satisfy the truncated set (|4.9p - (|4.12p 
for any postive k(9) so long as the twist vanishes, and for any k{9) > 3 if the twist is non zero. The 
value of A is irrelevant. 

Step 2. (Expansion with remainder terms.) 

As noted earlier, once we have the asymptotic solutions to the truncated system, we proceed to 
write the solutions to the Einstein equations as the sum of these asymptotic solutions plus remainder 
functions. Specifically, letting (v,w,y,(3) denote these remainders, we write 

U(6,t) = -(1 - k{0))lnt + <i>{9) + t £ v{6,t), (4.21) 
A(8, t) = ^(0) + t 2k{e) {B{9) + w{6, t)), (4.22) 

V(8,t) = \(l-k(6)) 2 \nt + \(9)+fy(9,t), (4.23) 
a(9,t) =j(9)+t E [3(9,t). (4.24) 

Step 3. (Fuchsian system.) 

We now show that if U,A,rj, and a satisfy the Einstein equations, then the remainder terms 
v,w,y, and (3 satisfy a system of the form 

(D + M{9))v = t £ ,f{t, 9, v, d e v). (4.25) 

If we can further show that, for analytic choices of k(8), (j>(9), ip(8), B(9), X(9) and 7(0) , we have (i) 
the nxn matrix M(9) satisfies the condition that, for any constant a S (0, 1), one has | \o M \ \ < C; and 
(ii) the function / is continuous in time, analytic in all other arguments, and extends continuously 
to t = 0, then we may apply Theorem 1 in [9] to prove that there exists a unique solution that is 
continuous in time, analytic in space, and vanishes as t — > 0. This would complete the verification 
of AVTD behavior. 
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Substituting the functions (|4.21[) - (|4.24p into the Einstein equations (|4.2|) - (|4.5p . and writing 
the result as a first order system by setting 

v = {vi,...,v & ) = (v, Dv, t E v e , w, Dw, t E w g ,/3, y), (4.26) 

we find that the PDE system satisfied by v takes the form 

D Vl =v 2 , (4.27) 

Dv 2 + 2ev 2 + e 2 Vi = t 2_£ (7 + t e v 7 )(-^k ee hit + <j) ee + dgv 3 ) 

+ (7 + fv 7 )[2k - e2t £ Vl - 2fv 2 )]t^ k2+3) - £ e 2( - x -' t ' }+2t€iv8 - Vl) A 

- il±^l[l(l — k) + Et £ Vl + t ^ 2 ]fl/2(fe-3)(fc+l)- £e 2(A+2t^ 8 ) K 2 

t 2-e 1 (4.28) 

H 2~ (79 + t £ dgv 7 )(--ke hit + 4>o + V3) 

f 2k-e 



2-2k-e 



2 

Dv 3 = t e do(ev 1 +v 2 ), (4.29) 
D«4 = v 5 , (4.30) 
Dug + 2fct- 5 = t 2 ~ 2fc (7 + ^7)^00 + t 2k {B + v 4 )((2k e hit) 2 + 2k gg \nt) 
+ At 2k k e \nt(B e + d 8 v A )] + t 2 {B 90 + w ee ){l + t £ v 7 ) 
- 2(7 + t e v 7 ){2k{B + Vi) + V5 ) t h(k 2 +3) e 2(\-4>)+2t'(vs~v 1 ) A 



(7 + t E V 7 ) t i/ 2 (k-3)(k+l) c 2(\+2t'vs) 
t 2-2k 



[2k(B + Vi ) + v 5 }K 2 



+ "-^-(W + t e 9ev 7 )[ipe + t 2k 2k e lnt(B + v 4 ) + t 2k (B e + d 6 v 4 )} 

+ 4t 2 ~ 2fe ( 7 + t £ v 7 )(ip g + t 2k 2k g \nt(B + v 4 ) + t 2k (B e + d e v 4 ))(-^kg lnt + <j) B + v 3 ) 

- 8kt £ (B + v 4 )(evi + v 2 ) - U £ v b {evx + v 5 ), 

(4.31) 

Dv 6 ^t'deievi + vs), (4.32) 
(D + e)v 7 = -4A( 7 + ^^K^-e^A-^t^-^) 

- (7 + f V7 f K 2 t l/2(k-3)(k+l)-e e 2X+2t'v 8 _ (4.33) 

(D + e)v s = t E {£ Vl + v 2 ) 2 + e(l - k)v 4 + (1 - k)v 2 
+ t 2 - E (j + t e v 7 ){--k 6 lnt + (/> e + v 3 ) 2 

+ L_ e ^+ 4 *^ (2k(B + v A ) + v 5 ) 2 



2-2k-e 



. e 4^f4f 0I (7 + t E v 7 ){ijg + t 2k 2k e \nt(B + v 4 ) + t 2k (B e + d e v 4 )) 2 



4 

+ A(7 + t € v 7 )t^ k2+3) ~ £ e 2{x ~' t '' >+2tS{v8 ' Vl) 

K 2 (^j + t e V 7 ) t i/ 2 (k- 3 )(k+l)-s c 2\+2t c v 8 

4 

(4.34) 

As written, the system (|4.27j) - (|4.34[) is of the Fuchsian form (|4.25p . However the right hand side 
does not satisfy the necessary conditions on /; that is, / does not extend continuously to t = 0. The 
equations for v 2 ,Vs, and «s have terms that would require k < 1. As stated earlier, terms involving 
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the twist constant K require k(9) > 3. All the terms requiring k < 1 involve either tjjg or ipgg, and 
can thus be eliminated by choosing ip(8) — constant in (|4.22[) . Thus if we assume -0 = constant and 
k(9) > 3, then / extends continuously to t = 0. The function / is also seen to be analytic in the 
other arguments, provided the functions k(9), <fi(9), and B{9) are analytic. The matrix M(9) is 

/ -1 00000 0\ 



M 



e 2 2e 

000000 

00 -1 000 

00 2/cOOO 

000000 

OOOOeO 



(4.35) 



\e(k-l) (fc-l) 00 ej 

That M(9) satisfies the bound | |ct m || < C for < a < 1, can been seen by direct computation. 

The argument just completed shows that for every choice of the analytic functions k(9), (f)(9), ip(9), 
B(9), X(9) and 7(6*) satisfying the conditions k(9) > 3 and ip(9) = constant, we have a T 2 -symmetric 
solution of the Einstein equations, either with A > or with A = 0, which has AVTD behavior. 
Comparing with the earlier work [9] on vacuum polarized T 2 -symmetric, we see that we have a wider 
class of solutions with AVTD behavior, since here we allow B(9) to be any analytic function, while 
in [9], we set it to zero. This result is consistent with the expectation that generic T 2 -symmetric 
spacetimes do not have AVTD behavior, but rather have the more complicated Mixmaster behavior. 
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